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Abstract. We study the appearance of a novel phenomenon for linearly coupled identical 
bursters: synchronized bursts where there are changes of spike synchrony within each burst. The 
examples we study are for normal form elliptic bursters where there is a periodic slow passage through 
a Bautin (codimension two degenerate Andronov-Hopf) bifurcation. This burster has a subcritical 
Andronov-Hopf bifurcation at the onset of repetitive spiking while end of burst occurs via a fold limit 
cycle bifurcation. We study synchronization behavior of two and three Bautin-type elliptic bursters 
for a linear direct coupling scheme. Burst synchronization is known to be prevalent behavior among 
such coupled bursters, while spike synchronization is more dependent on the details of the coupling. 

We note that higher order terms in the normal form that do not affect the behavior of a single 
burster can be responsible for changes in synchrony pattern; more precisely, we find within-burst 
synchrony changes associated with a turning point in the spiking frequency. 



1. Introduction. Elliptic bursting in a neuronal system is a recurrent alterna- 
tion between active phases (large amplitude oscillations) and quiescent phases (small 
amplitude oscillations). This kind of rhythmic pattern can be found in rodent trigem- 
inal neurons 16J, thalamic relay and reticularis neurons [5, 6J, the primary afferent 
neurons in the brain stem circuits [17j , and neurons in many other areas of the brain. 
It is clearly of interest for neuronal population information encoding and transmis- 
sion where several bursters fire within a population. Patterns of synchrony of elliptic 
bursters are may also be helpful in understanding firing patterns in more general types 
of burster H EE El US ■ 

In a previous study of the synchronization of elliptic bursters, Izhikevich exam- 
ined a pair of coupled "normal form" elliptic bursters [12] characterized by slow pas- 
sage through a Bautin (codimension two Andronov-Hopf) bifurcation. In that study 
burst (slow activity pattern) synchronization between the bursters was found to be 
easily achievable, whereas spike (fast activity pattern) synchronization was harder to 
achieve. Other studies include [7] who have examined nonlinearly coupled Bautin bi- 
furcations though not in a bursting setting and [ill El HI] wn0 have looked at various 
aspects of burst and spike synchronization for a variety of coupled burster models. 

In this article we study spike synchronization for coupled Bautin-type elliptic 
bursters with more complicated phase (spiking) dynamics. It transpires that higher 
order terms that are not important in the normal form of a single burster can be re- 
sponsible for nontrivial phase dynamics in coupled bursters. In particular, we observe 
and explain coexistence of and transitions between in-phase and anti-phase spiking 
within a single burst for two and more coupled bursters. This sheds light onto possible 
dynamical patterns of spike synchronization for coupled bursters in neuronal systems. 

We discuss a normal form for coupled Bautin-type elliptic bursters and focus 
on burst and spike synchronization in a system of n identical coupled bursters with 
zj e C, Uj e R and j = !,■■■ , n given by 
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(1.1) 



where cj, a, 77 £ R and 5, (7 € C are fixed parameters, i^j represents coupling. We set 

2 

B = 2 + zC = 2 + i^, C7=-l + *7=-l-i|. (1.2) 
We assume that the coupling term is 

n 

Kj = («i +i/c 2 )^CjfeZfe, (1.3) 

fe=X 

where ki,K2 € R are constant coupling parameters and Cjfe a constant connectivity 
matrix. For convenience here we take Cjk — 1 for j ^ k, Cjj = 0; i.e. all-to-all 
coupling. Biologically, although there are no rigorous reductions of specific bursters 
to this model, one can think of z = x + iy as a fast variable x that is analogous 
membrane voltage, y that is analogous to the fast current, and a slow variable u 
analogous to a slow adaptation current for a neuronal burster. 

The article is structured as follows: in Section [5] we discuss the individual burster 
behavior for the model equations p.l[) . In Section [3] we discuss two such bursters 
(n = 2), showing within-burst synchrony changes. These are analysed using a model 
system with assumed full burst synchrony in Section |4] and slow-fast dynamics |20j 
to reduce to an equation for within-burst phase difference. Bifurcation analyses of 
this equation help one understand the observed dynamics of the full model. Section [5] 
shows that we can observe similar dynamics in the three coupled elliptic burster system 
n = 3 and conclude with a discussion of some dynamical and biological implications 
in Section [6] 

2. The model for coupled Bautin bursters. Bursting is a multiple time scale 
phenomenon. In bursting, the fast dynamics of repetitive spiking is modulated by a 
slow dynamics of recurrent alternation between active and quiescent states. As ex- 
plained in [14] one may obtain bursting from a variety of dynamical mechanisms; here 
we focus on elliptic bursters (jl.ip with bursting behaviour by slow passage through a 
Bautin bifurcation; we briefly review the single burster dynamics. 

2.1. Normal form for Bautin bifurcation. Suppose we have a Bautin bifur- 
cation, namely a codimension two Andronov-Hopf bifurcation where the criticality 
changes on varying an additional parameter. Then there is a normal form that is lo- 
cally topologically equivalent to the bifurcation, and this normal form may be written 
[15] for z = x + iy 6 C as 

z = Az + Bz\z\ 2 + Cz\z\ 4 + 0{\zf), (2.1) 

where A, B, and C are complex coefficients. If we write A = A r + iAi, B = B r + iBi, 
and C — C r + iCi one can verify that an Andronov-Hopf bifurcation occurs as A r 
passes through and a change of criticality occurs where B r also passes through zero. 
The fourth order term is needed to determine the criticality at the degenerate point 
B r = 0. We write B and C as in (|1.2[) . It can be shown that Q, 7 and 0(|z| 6 ) terms do 
not affect the local branching dynamics of the system (|2.1[) . We will however argue 
that Q and 7 will influence the synchrony for two or more coupled elliptic bursters. 

From (|2.ip we obtain bursting dynamics [TTJ [121 [H] by coupling the system to 
a slow variable u € R that is the Andronov-Hopf parameter for the Bautin normal 
form, such that for z small u increases, while for z large u decreases: 

z = (it + iu)z + (2 + iC)z\z\ 2 + (-1 + i7)z|z| 4 
u = rj{a — \z\ 2 ). 

2 




Note that 77 <C 1 is the ratio of the fast to slow time scale. The system (|2.2p exhibits 
bursting for < a < 1 while tonic spiking sets in for a > 1. 
In polar form, z = re %e , (|2.2[) becomes 

r = ur + 2r 3 - r 5 ") 

(9 = w + Cr 2 +7r 4 i (2.3) 
m = r?(a — r 2 ). J 

In these coordinates it is clear that the fast subsystem undergoes an Andronov-Hopf 
bifurcation at u = and a limit cycle fold bifurcation (a saddle-node of limit cycles) 
at u = — 1. At the saddle- node bifurcation of limit cycles, stable and unstable limit 
cycles coalesce. A bifurcation sketch for the system (|2 .2[) is shown in the Figure [2~T1 
It is clear from this figure that the onset of periodic firing starts at a subcritical 
Andronov-Hopf bifurcation at u — with the emergence of a limit cycle. Likewise 
the steady state is reached via a saddle- node bifurcation of limit cycles at u = —1, 
where the stable limit cycle (solid line) meets the unstable limit cycle (dashed line) 
and eventually cancel each other at u = — 1. 

Note that during bursts, the limit cycles are non-isochronous; namely the inter- 
spike frequency 9 depends on r; there is a change in frequency of fast oscillation during 
the bursts. As this non-isochronicity does not affect the r or u dynamics, and hence 
the branching behaviour, it is not important for single bursters. The phase dynamics 
in (|2-3[) depends on amplitude r: 

= fi(r) = uj + (r 2 + 7 r 4 . (2.4) 




Fig. 2.1. Schematic bifurcation diagram for z in the system 12. 26 on varying u. SFP 
denotes the stable fixed point, UFP unstable fixed point, SPO stable periodic orbit and UPO 
unstable periodic orbit. It is clearly seen that at u = 0, the system undergoes subcritical 
Andronov-Hopf bifurcation, while saddle node bifurcation of limit cycles occur at U = — 1. 

The non-trivial periodic orbits of the system (|2.3j) are (ro,ito) = (\ // a 1 a 2 — 2a) 
for 77 = 0. Non-trivial periodic orbits, 7-^0 correspond to periodic orbits of (|2.2D 
with periodic spiking. The dynamics of the (|2.2p is summarized in figure 12.21 for 
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Fig. 2.2. Dynamics of a single compartment Bautin Burster governed by \2. 20 and $27, 
In panel (a), the timeseries of Re(z) is shown with solid line and the corresponding slow 
variable, u, with dashed line. The parameters for the simulation are u = 3, a — 0.8, n — 0.1, 
a = 4, and r m = 1.35. The arrows in (d) indicates the direction of change of the slow 
variable, u. 



parameters u> — 3, r\ — 0.1, a = 0.8, a = 2, (3 = —1, f = 0, and 7 = 0. Observe 
the slow passage effect [3] apparent from Figure 1^1 although the stability calculation 
shows the Andronov-Hopf bifurcation occurs at u = 0, but simulation shows a delayed 
bifurcation [3]. 

Note that we use parameters in (|1.2p such that 

C = — , 7 = -^, ( 2 -5) 

so that 

^ = ar(rl-ri). (2.6) 

For the parameters (o-,r m ) it is clear that there is a turning point at r — r m . The 
parameter a can be interpreted as the magnitude of non-isochronicity for the phase 
dynamics while r m is a turning point for f2 on changing r. 

2.2. Coupled elliptic bursters. We consider direct linear coupling of the nor- 
mal form system (|2.2I) via the fast variables z to give a coupled system of the form 
(jl.ll 11.3]) with coupling parameters Ki and K2- The coefficients Cjk for the coupling 
term of (|1.1[) are the connectivity matrix; here we assume all-to-all coupling, namely 



Cjk 



1 ttj^k 
otherwise 




Fig. 3.1. Within-burst synchrony change from stable inphase to stable antip hase st ates . 
This pattern repeats during each burst. This result is obtained from simulation of f 1. ij) , { 1.36 
and S1.3\) for n = 2, and parameters ki = 0.001, K2 = 0.2, cr = 3, r\ = 0.05, r m = 1.35, 
lo — 0.01. Noise of amplitude 10~ 5 was added to the fast subsystem. In this figure, the two 
coupled bursters are burst synchronized and the spikes become inphase at the beginning of 
the burst, but this synchrony pattern changes to antiphase near the middle of the burst. The 
inset in the topmost panel shows the region of the transition. Note that the initial transient 
and sudden change in the synchrony pattern along the burst profile are observable from d\2, 
where d\2 = indicates inphase synchronization. 



This form of coupling is analogous to the electrical (gap junction) coupling be- 
tween synapses with phase shift expressed by the argument of k\ + in2- Positive 
Ki corresponds to excitatory coupling, while negative k\ corresponds to inhibitory 
coupling. 

3. Burst and spike synchronization for two coupled bursters. We nu- 
merically investigate the dynamics of a pair of coupled elliptic bursters governed by 
the system (11.11) . Burst synchronization between the cells can be easily achieved for 
a wide range of parameter values with this system. In case of = and Ki > 
(excitatory coupling), this generally generates inphase bursts, while antiphase bursts 
result from inhibitory coupling^ 

There is a "within-burst synchrony change" observable within figure 13.11 The 
top panel shows X\ and x^. All transients were allowed to decay and the displayed 

1- We write the system l|l.lj) using z\ = x\ + iy± and Z2 = X2 + iy2 for the purposes of numerical 
simulation. All the simulations were done with the interactive package XPPAUT [8 t . For integrations, 
the built-in adaptive Runge-Kutta integrator was used, and results were checked using the adaptive 
Dormand-Prince integrator. 
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Fig. 3.2. Within-burst synchrony change from stable ant iphase to stable mphase states. 
The governing system and the parameters are as in fiaure [371\ except k,2 = —0.2. The inset in 
the topmost panel shows the transition in detail. In the last panel, the bump in di2 signifies 
the antiphase synchronization of the spikes within the synchronized burst. The corresponding 
slowly changing current variables, ui and U2, are shown in the middle panel with solid and 
dotted lines, respectively. 

pattern is repeated within each burst. A detail of the middle of the burst is shown in 
the top- right inset. The corresponding slow variables of the system, u\ and u 2 , are 
shown in the middle panel. The distinguishing solid and dashed traces correspond to 
the activity patterns of the two cells, respectively. The bottom shows the Euclidean 
distance 

di2 = y/ (xi - x 2 ) 2 + (yi - y 2 ) 2 + ("i - u 2 ) 2 , 

between the two systems to show the presence (di2 = 0) or absence {d\ 2 > 0) of 
synchrony. The values of the parameters used in the simulation are n,\ = 0.001, 
k 2 = 0.2, a = 3, r) = 0.05, r m = 1.35, u> — 0.01, so the two cells allowed to have 
same frequency. Wiener noise of amplitude 10 -5 was added to the system so that the 
system does not stick in any unstable state. 

The spikes are inphase at the beginning of the burst, but change to antiphase 
about the middle of the burst. The inset shows the region of this transition. This 
transition region may be shifted along the burst profile on changing r m . Larger values 
of r m shift this transition leftward along the burst profile, and vice versa. This sudden 
change in the synchrony pattern along the burst profile is also captured by di 2 . 

We present another example of within-burst synchrony change for different pa- 
rameter values in figure 13.21 where spikes of the two coupled cells start antiphase 
and change to inphase during the burst. The parameters for this are n\ = 0.001, 
k 2 = —0.2, r] = 0.05, cr = 3, and r m = 1.35. As before, a low amplitude noise of order 



10~ 5 was added to fast variables. The inset in the first panel shows the region of the 
transition. In the last panel, d\2 indicates that the burst is initially antiphase, and 
as it returns to d\2 = there is a transition to inphase synchronization of the within- 
burst spikes of the two cells. The corresponding slowly changing current variables, U\ 
and «2, are shown in the middle panel. The overlapped solid and dashed lines imply 
the inphase burst synchronization of the coupled system. 

4. A burst-synchronized constrained model. Although it is possible to find 
within-burst synchrony changes within , we are not able to explain their existence 
analytically from the model. To overcome this, we reduce the coupled system to a 
constrained problem where we assume burst synchronization, followed by a slow-fast 
decomposition. Using this we can explain how non-isochronicity and linear coupling 
can lead to within-burst synchrony changes. 

4.1. Two coupled bursters in polar coordinates. Writing (jl.ip with (|1.2p 



and (|1.3p in polar coordinates Zj = 



gives for n — 2 the system 



ri = u\rx + 2r\ - r\ + r 2 (/«i cos(6* 2 - 0±) - k 2 sin(0 2 - #i)) 

0\ = u + -ar 2 m r{ - -arf + sin(6» 2 - 9 X ) + k 2 cos(6> 2 - 6> x )) 

ui = r)(a - r\) 

r 2 = u 2 r 2 + 2r 2 - r 2 + n(«;i cos(6 l i - 6 2 ) ~ k 2 sin(6»i - 6 2 )) 



—i 

r 2 



u 2 = r)(a - r\). 



2 ) + k 2 cos(0i-0 2 )) 



(4.1) 



We constrain the system to exact burst synchronization by setting: 



u(t) = Ul(t) = u 2 (t) 



u = ui=u 2 = rj[a- -{\zi\ +\z 2 \) 



(4.2) 



Thus, the system (I4.1|) may be written with the constraint (|4.2p and considering 
-61 — 62 as 



UT2 
1 



2r\ 
2rl 



„5 



K\T2 COS ( 
KlTl COS ( 
1 



K 2 r 2 sm( 
K 2 ri sine 



= o°" r m( r i - r 2) - ~M r i - r 2) 



2 



r\ + r 2 



4 

sin <p — k 2 



rir 2 



COS ( 



(4.3) 



As we are interested in synchrony changes, we define longitudinal and transverse 
coordinates 



n = (r 1 +r 2 )/2 

n = (n - r 2 )/2. 
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(4.4) 



The system (I4.1[) reduces to the four dimensional system 



ri = uri 

+ K\Ti COS 

f t = ur t 

— K\Tt COS 
„2 



+ 2rf + 6nrf - rf - lOrf r\ - hnr\ * 



6rfr t 



t - 1 i 
K-2Xt sin cj) 

+ 2r\ - 5rfn - lOrfrf - rf 
K 2 ri sin 4> 

2 i „2- 



= 2ar m rir t - 2arir t (r l + r t ) 
-2^j + lysine 



4k, 



■ cos< 



it = 77(0 - (rf +r t 2 )). 



(4.5) 



Here, (ri,rt,4>) govern the fast dynamics, and u governs the slow dynamics. The 
system (|4. 5[) is a reduced four-dimensional realization of the full system (|4.1|) for a 
pair of coupled elliptic bursters. 

4.2. Stability analysis of the burst constrained system. In this section, we 
carry out a linear stability analysis of the fast sub-system of ()4.5|) about inphase and 
antiphase states with r t — and r; = r, which means both cells are burst synchronized 
with r\ = T2 = r. In the analysis, we assume the slow variable u is a constant of 
the system by setting the time scale ratio, 77, as a singularly perturbed parameter, 
i.e., 77 = 0. The dynamics, as a result, is only governed by the fast spiking activity. 
For r t = and 77 = we write r/ = r as the stable nontrivial solution of (|4.5|) in the 
appropriate subspace 



fi = (u + 2k\ cos 4>)ri + 2rf — rf 



(4.6) 



corresponding to bursting behaviour. Note that for small \ki |, this will have a solution 
close to the single burster case. 

If we consider the fast system of (|4.5[) with u between -1 and +1, then one can 
verify the existence of two solutions 

• Inphase where rt = 4> = 0, r; = r, 

• Antiphase where rj = 0, </> = 7r, 77 = r, 
where 7' is the solution of 



u = r 4 — 2r 2 — 2ki cos < 



(4.7) 



The Jacobian for the fast system at the inphase solution is block diagonal with a 
single real eigenvalue and a block 



u + 6r 2 — 5r 4 — K\ K^r 
2ar 2 m r - 2ar 3 - ^ -2k x 



(4.8) 



Likewise, the Jacobian for the fast system at the antiphase solution is also block 
diagonal with a single real eigenvalue and a block 



Ja 



u + 6r 2 — 5r 4 + K\ — K 2 r 
2ar 2 m r - 2ar 3 + ^ 2k\ 



(4.9) 



Note that the off-diagonal entries of the Jacobian matrices (|4.8j) and (J479J) depend on 
the imaginary part of the coupling coefficient, K2, and other system parameters. The 
real eigenvalues can be assumed negative because of stability of the solution of (14. 6|) . 



s 



The eigenvalues of the equation (|4.8p can be determined by examining the trace 
of the matrix (|4.8|) 

tr(J ln ) = u + 6r 2 - 5r 4 - 3/ti (4.10) 

and the determinant 

det(J in ) = -2(u + 6r 2 - 5r 4 )«i - 2err 2 (r 2 a - r 2 ) K2 + 2k\ + 4k 2 .. (4.11) 

Many interesting insights into the inphase dynamics of (|1.2|) and ()1.3|) may 

be extracted from (|4.10p and (|4.11|1 . Similarly, we can understand their antiphase 
counterparts from (14. 9p by examining 

tr(4»»«) = w + 6r 2 -5r 4 + 3Ki (4.12) 

and 

det{J anti ) = 2(u + 6r 2 - 5r 4 ) Kl + 2crr 2 (r 2 n - r 2 )n 2 + 2n\ + \k\. (4.13) 

For simplicity, we consider the special case of weak coupling where K\ — 0, and 
\k 2 \ -C 1. In such a case, it may easily be seen that both ir(Jj„) and tr(J ant i) 
in (|4.10p and (I4.12p . respectively, are negative, as stability of the periodic solution 
of (14. 6p means that u + 6r 2 - 5r 4 < 0. So, from (|4~TU|) and ^T2J, ir(J in ) < 0, 
and tr(J an ti) < 0. For this weak coupling, it is also evident that (tr( Jintanti))) > 
idet(Ji n t an ti\). Hence, the system will have a stable node for det(J in ( anti j) > and a 
saddle for det{J in(anti) ) < 0. 

To explain the within-burst synchrony change observed in figures 13.11 and 13 . 21 we 
write (|4.1ip and (|4.13p to first order in K2, with ki = and small k 2 > 0, as 

det(J in ) = -2crr 2 {r 2 n - r 2 )n 2 + 0{k 2 ) (4.14) 

and 

det(J anti ) = 2ar 2 {r 2 m - t 2 )k 2 + 0(k 2 2 ). (4.15) 

From equation (|4.14p . if r > r m + 0(k 2 ), then det(Ji n ) > 0. Together with the 
condition tr(Ji n ) < 0, this implies that the inphase solution is stable, whereas (|4. 1 5[) 
implies that the antiphase solution is unstable for r < r m + 0(k 2 ). We may derive 
approximate expressions for r(— r{) where the bifurcations take place. We denote the 
bifurcation value for amplitude of the inphase solution by r,- n , and the amplitude of 
the antiphase solution by r an u- Note that Ti n may be obtained by equating det(Ji n ) 
to zero in the equation (|4.1ip with K\ = giving 

det( J m ) = -K 2 {2ar 2 n - 2ar 2 m r 2 m - 4k 2 ) = 0. (4.16) 

Now solving JUTS}, 

r m =r m {l-^ r ) + 0{K 2 2 ). (4.17) 

Likewise, from equation (|4.13p . the bifurcation point, r anti , may be obtained as 
r an u = r rn (l + ^-) + 0(k 2 ). (4.18) 
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Fig. 4.1. Bifurcation diagram of 4> against u for the burst- synchronized contrained system 
f^.3| ), where u is a parameter that slowly decreases during each burst. The parameters are 
1.35, and different ki and K2 as indicated in the panels from (a) to 
(d). The solid lines represent stable solutions, while the unstable solutions are shown with 
dash-dotted lines. The arrow, running from right to left, shows the direction of the change 
of u within a single burst. 



Note that r depends on u via (|4.7|) . So, the corresponding bifurcation points for 
in(anti)phase oscillations can be derived from (14. 171 14. 18[) and (|4.7[) as 

u in = r 2 m (r 2 m - 2) + - r 2 m ) + 0(4), (4.19) 

and 

u anti = r 2 m (r 2 m - 2) - i^-(l - r 2 n ) + 0{k%). (4.20) 

4.3. Synchrony bifurcations of the fast system. We now consider numerical 
bifurcation analyses of the fast system (|4.3|) by taking r\ as the singular perturbation 
parameter to take the fast system through single bursts and to compare with the 
asymptotic results found for k\ = 0. 

To begin with we present in figure |4~T1 bifurcations of solutions of (|4.3|) projected 
onto the phase difference, 0, as u is varied. The solid line represents the stable so- 
lutions, while the unstable solutions are shown with dash-dotted lines. The arrow, 
running from right to left, shows the direction of the change of u. What figure [4~TT al 
shows is a burst that begins with stable inphase solution, and till almost half way 
through the burst, the inphase solution remains stable and then the antiphase so- 
lutions gain stability. The coupling coefficients in this results arc k\ — 0.001 and 
K2 = 0.2. The other parameter values are a = 3, u> = 3, and r m = 1.35. This be- 
haviour agrees with the simulation result shown in the figure ROI obtained from similar 
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Fig. 4.2. Two parameter bifurcation diagram of a against u for the fast subsystem of 
[4-3fy . The other parameters are fixed at «i = 0.001, K2 = 0.2, and r m — 1.35. There are 
stable inphase oscillations in region b,c, and stable antiphase oscillations in region a,b. 



k. =0.001 
K,=0.2 




-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

u 

Fig. 4.3. Two parameter bifurcation diagram of r m against u, for system and parameters 
as in figure \4-2\ and a = 3. There are stable inphase oscillations in the region b,c, and stable 
antiphase oscillations in the region a,b. 

system parameters. Similarly, figure 14.1( c) explains what is found in the simulation 
in figure EOl Here, the spikes in the burst start off in stable antiphase and changes to 
stable inphase. Figure I4~l7 b') and (d) show the results with k,\ — —0.001 but different 
K2- One interesting observation is the presence of the bistable region around the mid- 
dle of the burst separating the stable inphase and antiphase solutions. This region 
occurs near the transition point (r m = 1.35) along the burst profile as predicted in 
the analysis in the previous section. These bifurcations show the robust coexistence 
of the inphase and antiphase synchrony patterns of the within burst spikes for a range 
of u, and within-burst synchrony changes of the coupled bursting system 

Figure l4~2l shows continuation in the ua plane. The range of values of u is within 
burst activity period. As er increases, the bistable region is seen to get narrower, 
in agreement with (|4.17l I4.18p . Likewise, figure |4~31 is obtained from parameters: 
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Fig. 4.4. Bifurcation diagram of Ki against u for system and parameters as in fig- 
ure \4-l^ a,b) with «2 = 0.2. There are stable inphase oscillations in the region b,c, and 
antiphase in the region a,b. 
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Fig. 4. 5. Tw o parameter bifurcation diagram in Ki against u for system and parameters 
as in fjgure \4J\ c,d) with Ki varying and K2 = —0.2. There are stable antiphase oscillations 
in the region a,b, and inphase in the region b,c. 

K\ = 0.001, K2 = 0.2, and a = 3. What this result shows about the role of r m is 
that the position of the bistable region 'b' may be shifted along the burst profile by 
varying r m . 

The role of the coupling parameter ki is shown in figures [44] and [43] for two values 
of K2- The parameters in figure 14741 are K2 = 0.2, a = 3, and r m = 1.35. For these 
parameter values, the behavior is similar as in figure l4~TT a) and (b). It is interesting to 
note that within- burst synchrony changes appear even for weak inhibitory values (< 0) 
of K\. Moreover, stronger inhibitory values of K\ would mean only antiphase spike 
synchronization. Similarly, figure l4~5l demonstrates similar dynamics to figure [ITTT c) 
and (d). Figure [4751 has parameters as those in figure l4~4l except K2 — —0.2. The 
excursion of the bistable region 'b' above the dotted horizontal line indicates the 
appearance of within-burst synchrony changes for weak excitatory values, «i > 0. 
Stronger K\ results in inphase spike synchronization. 
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Table 4.1 

Comparison of the bifurcation points, rj n , r ant i, Ui n and u ant i, obtained from simulations of 
system \4-3ty and those from equations \4-lT\ \4- 18\ \4-19[ \4-®0V f or K i = 0, «2 = 0.2, <r = 3 and 
r m = 1.35. 





Tin 


Tanti 




Uanti 


From system (14. 3 J) (figure [474]) 
From equations (|4.17| 14.181 14.191 14. 20[) 


1.3210 
1.3229 


1.376 
1.3771 


-0.4433 
-0.4438 


-0.2027 
-0.2032 


Table 4.2 

Comparison of the bifurcation points as in table\4-l\ except «2 = —0.2. 




Tin 


T"anti 


Uin 


Uanti 


From system (|4.3|l ("figure 14. 5P 


1.376 
1.377 


1.321 
1.3229 


-0.2027 
-0.2032 


-0.4433 
-0.4438 


From equations (14.171 14.181 14.191 14.20(1 



Tables 14.11 and 14.21 show the comparison of the inphase and antiphase bifurcation 
points, ri n , Ui n and r ant i, u an u, for «i = and two values of K2 calculated from 
(|4.17l I4.19j) and (I4.18[ I4.20|) . respectively, with those from simulations of systems 
(|4.3p . Note that the bifurcation points obtained from the original system and those 
from the constrained system agree quite well. 

Bifurcation diagrams in figures 14.61 and 14.71 portray bifurcation diagram in UK2- 
space for two values K\. These figures show the role of K2 in spike synchronization. 
Figure 111)1 uses the parameters: K\ = 0.001, a = 3, r, m = 1.35. For negative and weak 
positive values of n-i (region c) only stable inphase solutions are present. Figure 14771 
shows the bifurcation diagram for k,\ = —0.001. As before a, b is the region of 
stable antiphase, and b,c the stable inphase solutions. It may be observed that for 
negative and weak positive values of K2 (region a), one would see only stable antiphase 
synchronization of spikes in the burst. 

To end this section we present in figure 14.81 an interesting result comparing the 
bifurcation point of the within burst synchrony change between the original system 
p. II 11.21 II. 3|) and the constrained system (|4.2p . Both systems have same coupling 
and system parameters. Note that the within burst synchrony change for the original 
system occurs at a more negative value of u than that of the constrained system. This 
difference may be attributed to the slow passage effect of the fast system within the 
burst dynamics [3]. We observed that the difference can be reduced by increasing 
noise amplitude added to the fast system. But eventually the large amplitude noise 
unsettles the stable solutions within the burst dynamics. 

5. Burst and spike synchronization for three bursters. We briefly demon- 
strate that within-burst synchrony changes are present in larger numbers of coupled 
bursters. In particular we look at three coupled Bautin-type elliptic bursters, i.e., 
(jl.ip . (jl.2p and (jl.3p with n = 3. A simulation of this system is shown in the fig- 
ure [54] with parameters to — 0.1, r m — 1.35, a = 5, n\ = —0.001, «2 = —0.2, and 
additive noise of amplitude 10~ 5 to the fast variables. 

The figure shows very similar behavior to the two-burster system with the dif- 
ference that the antiphase state is where the three bursters have a phase shift of 
relative to each other. The oscillations at the beginning of the burst are antiphase 
in this sense, and there is a transition to inphase during the burst. The inset in the 
top panel of the figure shows the transition in the spike synchrony pattern during the 
burst. The activity of the three different bursters is shown in solid, dashed and dash- 
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Fig. 4.6. Two parameter bifurcation diagram of K2 against u for system and parameters 
as in figure \^T\ a, b ) with n\ = 0.001 but ki varied. There are stable inphase oscillations in 
the region b,c, and antiphase in the region a,b. 
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-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 

I] 

Fig. 4.7. Bifurcation diagram as in figure \4-6\ but ki = —0.001 is varied. There are stable 
inphase oscillations in the region b,c, and antiphase in the region a,b. 

dotted lines. The middle panel shows evolution of the corresponding slow variables, 
Ui(t), U2(t) and U3(t). The third panel plots di2, e?i3 and c?23 that all must be zero 
for inphase synchronization, where 



dij = - Xj) 2 + {Hi ~ V]) 2 + (ui - Uj) 2 , 

with i, j = 1, 2, 3. 

6. Conclusion. We study the spiking dynamics of coupled elliptic bursters un- 
der direct linear coupling and find that within-burst synchrony changes are possible 
even for a simple normal form model, as long as terms that break isochronicity of 
the normal form are included. We observe that within-burst synchrony changes are 
stable and robust to changes in parameters. However, for identical bursters these 
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Fig. 4.8. Comparison of the bifurcation point of within burst synchrony change (a) for 
the two coupled system SI. 11 \1.2l \ 1.3(1 and (b) for the constrained system \4->ty - I n ( a )> 
the position of within burst synchrony change along u can be observed from drz (thick solid 
line) , where u = (ui + U2)/2, is defined as such to agree with constrained slow variable u in 
f4-5| ). The parameters for simulation in (a) are a = 3, r m = 1.35, ki = 0.001, K2 = 0.2, 
r) = 0.005, a — 0.8, and uj = 0.0003. A small amplitude noise of order 10 _J is added to the 
fast systems. Fast dynamics xi and X2 are shown in solid and dashed lines, respectively, (b) 
depicts the inphase bifurcation solution branch of the fjgure \^.l\ a) . The solid line represents 
the stable inphase solution, and the unstable inphase solution is shown by dashed line. Note 
the difference in the inphase bifurcation points, Ui n , and the loss of stability within burst. 



within-burst changes are only easy to observe in the presence of noise; this is because 
in the absence of noise the system may become stuck in unstable synchronized states. 

By reduction to fast-slow dynamics for the constrained burst-synchronised model 
we analyse the appearance of the within-burst synchrony change for two oscillators, 
and the influence of various system parameters. In particular we find that a turning 
point in the frequency can be associated with the observed within-burst synchrony 
changes, analogous to bifurcations observed in systems of coupled weakly dissipative 
oscillators [2J- Moreover, we can find the approximate location of the transition be- 
tween stable inphase and antiphase oscillations from bifurcation analysis of a reduced 
system. It will be a challenge to generalize this analysis to a case where the system 
is not burst-synchrony constrained. The examples we have illustrated in this paper 
are clearest for "long" bursts where there are many oscillations during which the syn- 
chrony changes. Similar effects are presumably also present in shorter bursts, but are 
harder to observe because the changes in synchrony must occur over a small number 
of spikes to be observable. 

For larger populations of oscillators we expect there can be not just transitions 
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Fig. 5.1. Within-burst synchrony change from stable antiphase to stable mphase states 
for three coupled Bautin-type elliptic bursters; see text for details. All three cells burst syn- 
chronously, but the fast spikes are antiphase at the beginning and inphase by the end of the 
burst. 

between inphase and antiphase during bursts, but also spontaneous changes in clus- 
tering, leading to robust but sensitive phase dynamics [HUB] and we believe this study 
gives some insight into the range of synchrony dynamics of coupled bursters in gen- 
eral. Better understanding of spike synchronization in more general coupled burster 
networks may lead to better understanding of potentially important new mechanisms 
for information processing and transmission by coupled neuronal bursters. This is 
discussed for example in [13j where it is suggested that information transmission may 
occur via resonance between burst frequency and subthreshold oscillations. 
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